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Abstract. We calculate the force operator for a charged particle in the field of an 
Aharonov-Bohm flux line. Formally this is the Lorentz force, with the magnetic field 
operator modified to include quantum corrections due to anomolous commutation 
relations. For stationary states, the magnitude of the force is proportional to the 
product of the wavenumber k with the amplitudes of the 'pinioned' components, the 
two angular momentum components whose azimuthal quantum numbers are closest to 
the flux parameter a. The direction of the force depends on the relative phase of the 
pinioned components. For paraxial beams, the transverse component of our expression 
gives an exact version of Shclankov's formula [Shelankov A 1998 Europhys. Lett. 43, 
623 - 8] , while the longitudinal component gives the force along the beam. 

Nonstationary states are treated by integrating the force operator in time to obtain 
the impulse operator. Expectation values of the impulse are calculated for two kinds of 
wavepackets. For slow wavepackets, which spread faster than they move, the impulse is 
inversely proportional to the distance from the flux line. For fast wavepackets, which 
spread only negligibly before their closest approach to the flux line, the impulse is 
proportional to the probability density transverse to the incident direction evaluated at 
the flux line. In this case, the transverse component of the impulse gives a wavepacket 
analogue of Shelankov's formula. The direction of the impulse for both kinds of 
wavepackets is flux dependent. 

We give two derivations of the force and impulse operators, the first a simple 
derivation based on formal arguments, and the second a rigorous calculation of 
wavepacket expectation values. We also show that the same expressions for the force 
and impulse are obtained if the flux line is enclosed in an impenetrable cylinder, or 
distributed uniformly over a flux cylinder, in the limit that the radius of the cylinder 
goes to zero. 
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1. Introduction 

There have been a number of investigations of the force exerted on a charged particle 
by an Aharonov-Bohm flux hne. Classically, of course, there is no force, so it, like the 
Aharonov-Bohm effect itself, is essentially quantum mechanical, vanishing as h ^ 0. 
Olariu & Popescu (1983, 1985) show that for certain localized wavepackets (these are the 
fast wavepackets described in Section below), the force, along with the momentum it 
imparts, is negligible unless the centre of the wavepacket hits the flux line. Nielson and 
Hedegard (1995) and Shelankov (2000) compute matrix elements of the force operator 
for stationary states of the same energy. Shelankov (1998) calculates the transverse force 
on a stationary incident beam of finite angular width using a paraxial approximation, 
a result we refer to as Shelankov 's formula. This use of the paraxial approximation has 
been justified by Berry (1999), who computes the asymptotic deflection of the beam 
from an exact representation. Peshkin (1981, 1989) computes the expectation value of 
the force when the flux line is enclosed in an impenetrable cylinder. Recent interest in 
this problem has been stimulated by the analogy with the lordanskii force (lordanskii 
1966) exerted on phonons by a vortex in a superfiuid (see, eg, Sonin 1975, Sonin 1997, 
Thouless et al 1997, Stone 1999), which has been the subject of some recent debate. 

In this paper we add to these investigations in several ways. First, we obtain 
exact expressions for the Lorentz force operator due to an Aharonov-Bohm flux line 
and its matrix elements between arbitrary stationary states; the restriction to on-shell 
matrix elements, between stationary states of the same energy, agrees with previous 
results. We also obtain an exact expression for the time integral of the force, the 
impulse operator, and compute its expectation values for various kinds of wavepackets. 
We show that the force operator can be simply derived, formally at least, from purely 
kinematic considerations, and give a mathematically rigorous demonstration to justify 
the results obtained in this way. 

The paper is organized as follows. In Section |^, we give a formal derivation of 
the (vector) Lorentz force operator due to an Aharonov-Bohm flux line. Pointing out 
that the nominal magnetic field operator, {ahc/e)S'^{r), is incompatible with gauge 
invariance, we show that that a modification of the canonical commutation relations 
restores gauge invariance and yields the Lorentz force. An explicit formula for the 
expectation value for stationary states follows directly. On-shell matrix elements 
between stationary states of the same energy are seen to coincide with previous results. 
The transverse component of the force is shown to coincide with Shelankov's formula in 
the paraxial limit (see also Shelankov (2000)). We then compute (Section ^ the impulse 
operator, the integral of the force operator in infinite forward and backward time, in 
the position representation. In Section ^ we compute the leading-order expectation 
value of the impulse for two kinds of wavepackets. For slow wavepackets (Section [4.1| ), 
which spread more quickly than they move, the impulse is inversely proportional to 
the initial distance between the wavepacket and the flux line, with magnitude and 
direction depending periodically on the flux parameter a. By treating the deflection of 
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a slow wavepacket as a classical scattering, we obtain an expression for the scattering 
cross-section which, surprisingly, coincides with the exact result. Fast wavepackets 
(Section |4.2| ) move more quickly than they spread, so much so that the fractional increase 
in their width as they pass the flux line is small. In this case, the impulse is proportional 
to the transverse probability density at the flux line, and its transverse component gives 
a wavepacket analogue of Shelankov's formula. 

In the Appendix, we give a rigorous derivation of the force and impulse expectation 
values for a class of well-behaved wavefunctions. In Section ^, we compute matrix 
elements of the force operator for two standard regularizations of the Aharonov-Bohm 
flux line, the first where the flux is enclosed in an impenentrable cylinder, and the second 
where it is distributed uniformly in a cylindrical tube. The Aharonov-Bohm results of 
Section |^ and Section |^ are recovered in the limit that the radius of the cylinder goes 
to zero. 



2. Force operator 

Consider a particle of charge e and mass M moving in the xy-plane in a magnetic 
field along z. Quantum mechanically, the particle is described by the Hamiltonian 
H = \MV ■ V, where MV = P — cA/c is the kinetic momentum and A{r) is the 
vector potential. The force, ie the rate of change of the kinetic momentum, is given by 
the appropriately symmetrized Lorentz force operator, 

F = i[MF, H] = ^^{VA{Bz)- {B z)AV), (1) 

where the magnetic field operator is defined by 

B = ^[MV^,MVy]. (2) 

Here $o = ^Tchc/e is the magnetic flux quantum. The usual commutation relations for 
position and canonical momentum lead to the usual relation between the magnetic field 
and the vector potential, namely B z = S/ A A. 

However, this relation is incorrect for the vector potential of an Aharonov-Bohm 
flux line (since we are restricting to the plane, we should perhaps say 'flux point', but 
we will follow conventional usage). For a flux line at the origin of strength a$0; and in 
a circularly symmetric gauge, the vector potential is given by 

Air)=a^J. (3) 

As is well known, physically meaningful quantities depend only on the fractional part, 
a, of the flux parameter, a; a unit shift in the flux parameter, a — >■ a -|- 1, is equivalent 
to the gauge transformation ip — > Uip, where 

{Uij){r,^)=e'U{rA)- (4) 
As a consequence, a physically observable operator 0{a) which depends on the flux 
parameter must transform under U according to 

UO{a)U^ = 0{a + l). (5) 
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For example, the kinetic momentum MV satisfies this relation, as is easily verified. As 
the magnetic field operator is expressed in terms of the commutator of components 
of kinetic momentum, it must also satisfy (^. However, V A A, which is given by 
a$o^^(^) ^1 clearly does not satisfy (|^); it is invariant under the gauge transformation 
U (like any operator which depends only on position), but is not periodic in a. It follows 
that substituting V A A for the magnetic field in (|I|) cannot give the correct expression 
for the Lorentz force operator. 

The problem is caused, of course, by the singularity at r = 0, and can be avoided by 
an explicit calculation of the time rate of change of the expectation value of the kinetic 
momentum for suitably chosen wavefunctions. This is done in the Appendix. However, 
we can obtain the same result more easily and more directly from formal arguments. If 
the usual canonical commutation relations lead to a magnetic field operator which does 
not transform correctly under gauge transformations, then the canonical commutation 
relations must be modified by the fiux line. In particular, we show below that, formally, 
the components of the canonical momentum, and Py, do not commute - equivalently, 
the partial derivatives and dy do not commute - in the presence of nonzero fiux. 

First, we note that the partial derivatives dx and dy certainly commute when applied 
to smooth wavefunctions. That is, 

[5.,5,]^ = (6) 

for '4){r) smooth. Applying the gauge transformation U to this relation m times, we get 

[U'^dJJ^'", W^dyU^""] {e'"">'^) = 0. (7) 
The partial derivatives transform according to 

W^dxU^" = dx- imdx^, W^dyU^" = dy- imdycj). (8) 
Substituting (|]) into (|^), and using the differential version of Stokes' theorem, 

[a,,9,]0=(VAV),0 = 27r52(r), (9) 

we obtain 

[dx, dy] (e''"'^V') = 27r(52(r)ime''"V = 27r(5'(r)a<^ (e''"^) , (10) 

where the second equality follows because d^ip vanishes at the origin for il){r) smooth. 
This implies, formally, the operator relation 

[dx,dy\=2T,5\r)d^, (11) 

or, equivalently, 

[Px,Py] = -2Txm\r)L, (12) 

where L = {h/i)d^ is the canonical angular momentum. 

Using the modified commutation relation ([T^ to evaluate the magnetic field 
operator (^, we get, instead of the classical relation B = a^o6'^{r), the result 

^ = 2^[^^'' + = "T^'^^^^' ^^^^ 
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where 

Az = r AMV = {L~ah)z (14) 

is the kinetic angular momentum. It is evident that A, and hence B, satisfy the 
transformation law (|^). 

Throughout this and the following sections, it will be convenient to represent vectors 
in the xy-plane as complex scalars. For example, W = x + Wy y will be represented 
by W = Wx + iWy. If W is a vector of hermitian operators (as opposed to real scalars), 
then W is the nonhermitian scalar operator whose hermitian and antihermitian parts 
are Wx and iWy respectively. In this way, the kinetic momentum MV is represented by 
the scalar operator 



MV 



h 



ids + a 



(15) 



It is useful to note that, in general, z f\W is represented by iW. 

We proceed to compute the expectation value of the Lorentz force. Substituting 
(p!3|) for the magnetic field and (|15]) for the kinetic momentum into the expression (|lD, 

we get 



^1 



e 

-1- 

c 



BV\ 



2i^ 



2ix 



6{r] 



la e 



ids + a 



ilj{r,(f))rdr, (16) 



where we have used 6'^{r) = 6{r)/{7[r). With ip resolved into its angular momentum 
components. 



im(f> 



(17) 



m=— oo 



the integrals in (0) are trivially evaluated (note that S{r) dr = i), with the result 



2^' 



m=— oo 



[m — a, 



J T- = 



Like any vector operator, JF couples only consecutive angular momentum components, m 
and m + 1, and, as one would expect, depends only on the behaviour of the wavefunction 
at the flux line. For well-behaved wavef unctions, the leading-order behaviour of ipmif') 
as r — > is given by 

V'm(r)~C^rl'"-°l. (19) 

Sufficient conditions for (|TUp are discussed in the Appendix. Here, we note that (^) 
ensures that the energy density, ip*ii^)iHip){r), is finite at r = 0. Substituting (|T9D into 
(H), we get 

\^\^) = -^rr- (m + 1 — a)(|m — a| — (m — a)) 



M 

m=—oo 

/-I* |m+l-Q:| + |m— a| — 1 



r=0 



(20) 
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The terms in the sum (^) vanish unless m = a, where a = a — a denotes the integer 
part of the flux parameter. Thus, only the 'pinioned' components of the wavefunction, 
ipa and ipa+i ^ those whose angular momentum quantum numbers are nearest the flux 
parameter a - contribute to the force expectation value (p^S]). From (pOD, 

(^|^|^) = ^«(l-a)C:+,a, (21) 
or, equivalent ly, 

4:7!- 

J" = ^a{l - a) \Ui) iU , (22) 
where the state |^m) corresponds to the singular wavefunction 



6{r) 



^"^(-''^) = ;^7^«'"^^ (23) 

so that {^m li') = Cm- It is readily verified that the force operator ( p2D transforms 
according to (|^) under the gauge transformation (^). A rigorous derivation of the 
expectation value ( pT]) for suitably chosen wavefunctions is given in the Appendix. 

In the preceding derivation, the force operator due to a flux line, like the force for 
a nonsingular potential, is derived from kinematics, specifically from the commutation 
relations. The derivation does not require the solution of the Schrodinger equation. 
Thus, it is straightforward to generalize to the case of more than one flux line (for which 
solutions of the Schrodinger equation are, in general, not available); the force operator 
is just a sum of contributions (^) centred around each flux line. 

However, to calculate the force on stationary states, or the time dependence of the 
force on nonstationary states, it is necessary to solve the Schrodinger equation. As is 
well known, eigenstates of the Aharonov-Bohm Hamiltonian with energy h^k'^ /2M and 
angular momentum mh are given by 

XkAr) = J\^-a\{kry'^t (24) 

where Jp{z) is a Bessel function. From the small- 2; behaviour, Ju{z) ~ {z/2Y /Tiy + 1), 
and the reflection formula, T{v)T{l — v) = tx / sinvrz/, we obtain from (p2D the matrix 
elements 

(Xp,nl-^IXfe,m) = °'^Vinia5m,Jn,a+l- (25) 

For = p, ie for stationary states with the same energy, we obtain the on-shell matrix 
elements 

{XKn\^\XKm) = —ksinna 5m,aSn,a+l, (26) 

in agreement with results of Nielson and Hedegard (1995) and Shelankov (2000). 

A general stationary state |\E') is a superposition of eigenstates with k fixed, and 
may be taken to be of the form 

CO 

l*)= E i-^y^""^bm\Xk,m) . (27) 
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($1 = __e-i™fcsin7ra6;;+i6a. (28) 



For bm = 1, \&(r) corresponds to a scattered plane wave incident from the right 
(Aharonov & Bohm 1959). From (pOf), we get the expectation value 

Shelankov (1998) has obtained an approximate formula for the transverse 
component of the force acting on a stationary beam of finite angular width. His analysis 
is carried out in a singular gauge, in which the vector potential vanishes everywhere 
except along the ?/-axis. A stationary beam incident from the right (in fact, Shelankov 
takes the beam incident from the left, but we revert to the convention of Aharonov and 
Bohm) is taken to be of the form e~^^^'il){x,y). Treating ip^^ as small compared to k'll)^ 
amounts to a paraxial approximation, in which the wave evolves freely in x (with x 
playing the role of time) for x 7^ 0, and is scattered by the vector potential at x = 0. 
The change Apj^ in the transverse kinematic momentum, {h/i) J°°^ilj*{x,y)'iljy{x,y) dy, 
is then calculated to be 

|2 



Apy = h Sin 27ra ^f'^^^^^ (29) 
j_^\i^in{y)\ dy 



V'in(z/) = ^^/ a{ky)e"'yy dky (30) 



00 



where 



/27Tk J-c 

is the incident wave at x = 0^, expressed here in terms of its transverse Fourier 
amplitudes a{ky). Multiplying Apy by the incident flux, which is given paraxially by 
{hk/M) |^m(y)|^ dy, gives Shelankov's formula for the transverse force, 

Ff) = ^fcsin27r5|^i,(0)^ (31) 

We now show that the y-component of the exact force expectation value, ie the 
imaginary part of (PE|), coincides with Shelankov's formula (^) in the paraxial regime. 
(Shelankov (2000) gives the same argument.) As discussed by Berry (1999), the state 
(^) can alternatively be viewed as a superposition of scattered waves incident from the 
directions (cos 6', — sin 6'), with amplitudes A{9) related to the coefficients b^ according 
to 

1= [ A(0)e^('^-")^ d0. (32) 
27r J-TT 



The paraxial approximation is valid for A{6) strongly peaked around = 0, with angular 
width w << 1. In this case. Berry (1999) has shown that A{6) ~ a{k6). From (pO|) 
and (|32D, it then follows that bm ~ '^in(("^ — a)/k) for |m — a\ « l/w, so that 

&:+l&a~|^in(0)p. 

3. Impulse operator 



For nonstationary wavepackets ipir), whose wavefunctions are not eigenfunctions of 
the Hamiltonian, the expectation value of the force does not itself have much physical 
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significance. It depends on the behaviour of the wavefunction near the flux hne, 
regardless of where the wavepacket is localized, and can oscillate rapidly as the 
wavepacket evolves. Of greater physical interest is the impulse imparted to the particle 
over the course of its evolution, either in the past or future. Let 

J^(t) = e^^*/^J^e-'^*/^ MV(t) = e'^*/^MVe-'^*/^ (33) 

denote the time-evolved force and kinetic momentum operators. Then MV{t) = J^{t). 
The forward (+) and backward (-) impulse operators are defined by 

= MV{±oo) - MV(0) = / J^(t) dt. (34) 

Jo 

Let us compute the kernal of the impulse operator in the position representation, 
I±{s,r) = (s|Xj- \r). From the completeness relation, 

^ $Z / \Xk,m) {Xk,m\kdk = 1, (35) 

we obtain 



2tt 

m=—oo 







1 p±oo 



m=— oo n=~oo 







X 



noo poo 

/ / \Xp,n) {Xp,n\ ^ \Xk,m) {Xk,m \r) exp(i^(p^ - k'^)t/2M) kdkpdp. (36) 
^0 Jo 



From the expression (^51) for the matrix elements {Xp,n\ J^\Xk,m), the only contribution 
to the double sum in (|36D is from the term n = a + 1, m = a. Substituting the 
eigenfunctions (p^, and letting {s,9) denote the polar coordinates of s, we get 



oo 



_ , , h"^ sirnia ... ^ , „ . ,x ; , 
J±(s, r) = ± ^^^^ exp(i(a + 1)6 - ia0) / dt 







X ^"exp (±i^t) Ji_5(p.)p2-"dp^"exp (ri^t^ Ukr)k'+^ dk. (37) 
The k- and p-integrals are of the form (Abramowitz & Stegun 1970) 



oo TV 



e- " Uhu)u'^^' du = ^^^e-^ . (3^ 
Substituting this result into (^), we get 
J±(s, r) = " sin ira exp(i(a + 1)6* — ia0 =F ivra)s"'^^"r" 

With the substitution t = 1/w, the remaining integral in (pQ]) is of the elementary form 

^ 'V ' ^'' ) = (0-^ ± i(H - .^))^ - 



with cr^ = 2;i/M. Substituting (H) into (H), we obtain 

2^±(s,^') = ^sinvra exp(i(a + 1)^ - ia0 ^ ivra)— — — — . (41) 
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One can verify that this expression transforms correctly, ie according to (^), under the 
gauge transformation (^). 

The denominator in (^Tj) can be alternatively expressed as 

(0+^i(,l_,2))2 = V(r^ - s')) ± ^rrS'ir^ " s% (42) 
Here P'{l/x), the derivative of the principal part, acts on functions f{x) according to 

POO 

f{x)P\l/x) dx = - /^dd(a:)/xda;, (43) 



where /odd(^) — i(/'(^) ~" /'(^^)) denotes the odd part of f'{x). For subsequent 
calculations, however, the integral representation ( ^0[ ) will be more convenient. 



4. Expectation values of impulse for wavepackets 

We parameterize wavepackets by their position R, width a and kinetic momentum hk. 
A convenient form is 

1 fr-R\ 

^{r) = -f exp(-ifc • r + ia0). (44) 



a \ a 



Here f{u) is a smooth normalized function localized at the origin with unit width and 
vanishing average (dimensionless) momentum, ie 



rfd'u = i,j J rfud'u = o, jjrfu'd'u = i, jjrvfd'u = o. (45) 

We assume that a << i?, so that the wavepacket ip{r) is localized far from the flux 
line. The phase factor exp{ia(j)) in (|4^) insures that hk is the kinetic, rather than the 
canonical, momentum of the wavepacket; its branch is chosen so that the phase factor 
is continuous over the region where ip{r) is appreciable. 

From (^) and (0), the expectation value of the impulse is given by 



2\h 

(J±)(i?,fc,a,«) = {tP\I±\tP) = — -e^^^sinvra 




f* I £ —\ f ( - ?^ e^Mr^^^ i(l-a)e+ia</, a 1-5 j2^2 ^Q) 



Since / has unit width, the integrand in (|46|) is appreciable only for |s — i?| ~ cr and 
\r — R\ cr. In this region, we can, to leading order in cr/i?, replace the phase factors 
exp(i(l — Q.)0) and exp(ia0) by exp(i(l — (5)$) and exp(ia$), respectively, where $ is 
the polar angle of R. Likewise, we can replace the factor r^s^""" by R. With the change 
of variables u = {r — R)/a and v = {s — R)/a and the integral representation (|40|), 
(^) becomes 

(X±) = — — e + smvra 

R u 



CO 



X / dw w 







/(n) exp ( icrfc ■ u =]= 2iw =F iwu J d u 



(47) 
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Thus, the direction of the impulse, arg is given by 

arg (J±) = $ + (i T a)vr. (48) 

For a = ^, the forward impulse is directed away from the flux line, and the backwards 
impulse towards the flux line. 

There are two parameter regimes where the expression (^7p has a simple asymptotic 
form, namely ka << 1, which corresponds to slow wavepackets, and ka » R/a, which 
corresponds to fast wavepackets. These cases are discussed separately below. 



4-1- Slow wavepackets. 



The condition ka << 1 implies that the wavepacket spreads (with velocity ~ h/Ma) 
more quickly than it moves (with velocity hk/M). Since / has unit width, the integrand 
in (^) is appreciable only for u of order one. In this case, for ka << 1, the phase factor 
exp{iak-u) is nearly equal to one. On the other hand, the phase factor exp{^2iwR-u/a) 
oscillates rapidly in this region, and hence renders the integral small, unless w is small, 
of order a/R. For w of order a/R, the phase factor exp{m'^w) is nearly equal to one for 
u of order one. To leading order in a/R and ka, ( ^Tf ) becomes 

2ih R 



a^ 



sm Tta 



dw w 



/(n)e^2™^-"/" d' 



u 



sm na 



w dw. 



where 



denotes the normalized Fourier transform of f{u). Letting 



Pie) 



hie) 



(49) 



(50) 



(51) 



denote the probability distribution for the direction, e, of the dimensionless momentum, 
^ = ^ e, we can write 

21^1 



(J^) = =^e'*^^™sin7rap(±^). 
R 



(52) 



Note that if f{u) is circularly symmetric, then p( e) is equal to l/27r. 

To leading order in a/R and ka, the impulse ( |52D is independent of the width 
and momentum of the wavepacket, and is of order h/R (ie, inversely proportional 
to the distance from the flux line). This is much smaller than the dispersion of the 
momentum, which is of order h/a. Therefore, to detect the impulse on slow wavepackets 
experimentally, one would have to perform a large number of measurements (on the 
order of (R/a)"^) of the asymptotic momentum on an ensemble of identically prepared 
systems. 
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Figure 1. At i = the wavepacket is centred ai Y y and moves in the — x direction. 
(3± are the angles of the incoming and outgoing velocities with respect to — i. 



By treating the motion of the centre of a slow wavepacket as a classical trajectory, 
we can derive an expression for the scattering cross-section a (6) using the classical 
formula, 



-{0) 



(53) 



Here h is the impact parameter, and 9 is the angular direction of the outgoing trajectory. 
Consider a slow wavepacket (0) centred on the |/-axis at F ^ at t = (thus, i? = |y| 
and $ = sgn(y)7r/2), moving in the —x direction with kinetic momentum hk. For 
simplicity, we take f[u) to be circularly symmetric, so that p = l/27r. Let /3_ denote 
the angle between the velocities at t = —oo and t = 0, and /3+ the angle between the 

= oo (see Fig 1). From (|5^) , these are given by 

Ke {ip\X±\ilj) — hk siniT a cos IT a — TckY 



velocities at t = and t 

cot /3± = : 



(54) 



Im {'il^\I± IV') sin vra 

(thus /3_ = /3+). Because the impulse is circularly symmetric, the angles /3± are 
unchanged if we rotate the entire system so that the incoming velocity, at t = — oo, 
is in the — x direction. In this case, the direction of the outgoing beam is given by 

^ = 7r + /3+ + /3_ = 7r + 2/3+. (55) 

To determine the impact parameter b, we appeal to classical angular momentum 
conservation, MV^b = sgnb MVoRq, where V- is the speed at t = —oo, and Ro and Vq 
are the distance and speed at the point of closest approach to the flux line. For the 
Aharonov-Bohm Hamiltonian (and, indeed, for any purely magnetic Hamiltonian) , the 
speed V = is a conserved quantity. Thus b = sgnbRo. We take the point of 

closest approach to occur at t = (when the velocity of the wavepacket is orthogonal 
to its position), so that b = Y. Then, from (|53|) , (|5^) and (^5|), 



aie) 



db 




d9 


-1 




-1 


de 




db 




dV 





2 sin^ (3. 



d(cot 
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27rfccos2^/2' ^ ^ 

Surprisingly, the expression (|5^) agrees with the exact result found by Aharonov 
& Bohm (1959). Of course, the preceding should not be regarded as a legitimate 
derivation of the scattering cross-section. Apart from certain ad hoc elements (eg, 
circularly symmetric wavepacket and the determination of the impact parameter), a 
proper derivation of the cross-section from time-dependent solutions of the Schrodinger 
equation requires wavepackets (unlike the slow ones used here) whose momentum is 
sharp. Still, it is interesting to ask whether or not this agreement is purely fortuitous. 



4-2. Fast wave packets 



A wavepacket initially at a distance R from the flux line with kinetic momentum hk 
reaches its point of closest approach to the flux line after a time r of order R/{hk/M). 
It spreads at a speed W of order h/Ma. Thus, at closest approach it will have spread 
a distance of order Wt ~ R/ka. For this to be much less than the width a, we require 
R << ka"^, which is just the condition for fast wavepackets. 
Let 

2 



J J f(u) exp (iak ■ u =]= 2\wR ■ uja ^ iwu'^) d" 



u 



(57) 



denote the u-integral which appears in the impulse expectation value (^TQ . Because / 
is of unit width, the integrand is appreciable only for u of order 1. For ka » R/cr, the 
phase factor exp[iak ■ u) in (^7|) is rapidly oscillating, and hence leads to a vanishingly 
small integral unless it is balanced by the phase factor exp(=pi2«wi2 ■ u/a). For such a 
balancing to take place, w must be large, of order ka'^/R. Therefore, the quadratic phase 
factor exp(=i=iw-u^) is rapidly oscillating, so ([57|) can be evaluated using the stationary 
phase approximation. To leading order in 1/w ~ R/lka"^), we obtain 



f 



To^k/2w - R 



a 



{51 



Substituting into (|i7|), we get 

R 



2in^e'*^^™ sinvra 



2ini?e^*^^^° sin 



R 



il){^r k) 



a 

2 

r 



dw 

w 



(59) 



where we have used (p^ to express the integral in terms of the wavefunction ip{r). 
Note that ip{r) behaves for small r as r" or r^~" (cf ([T9|)), so that the integral in (^) 
is convergent. 

In what follows, let us assume for concreteness that k is directed along - x, so that 
the wavepacket is moving to the left. We write R = X x + Y y. Unless the wavepacket 
is centred near the x-axis (specifically, unless \Y\ ~ a), ■ip{^rk) will be negligible over 
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the range of integration in (|59|). Thus, to leading order in cr/R, we may take i?e'* = X. 
Substituting this result into (0), we obtain the expression 

(J±) = 2i/ie^'™Xsin7ra / \ij{x,0)f —. (60) 



X 

for the impulse. 

Since the wavepacket is centred near X x, the x-integral in (pDf) is negligible unless 



X > in the forward (+) case (so that the wavepacket is moving towards the flux line), 
or unless X < in the backward (-) case (so that the wavepacket is moving away from 
the flux line). Assuming that ±X > 0, the main contribution to the integral comes 
from |x — X| ~ a, so that, to leading order in a/R, we can replace the factor 1/x by 
1/X in (0), and extend the lower limit of the x- integral to ^oo. Letting 

POO 

i^trans(y) = / dx\ij{x,y)f (61) 



denote the wavepacket 's probability density in y (the direction transverse to the incident 
velocity), we obtain, to leading order in a/R and R/lka"^), the expression 

(X±) = ±2ineT^™sin7r«0(±X)Ptrans(O) (62) 

for the impulse on fast wavepackets. Here G(x) is the unit step function. 

The impulse (p2) is independent of the wavenumber k. To leading order, it vanishes 



for wavepackets which miss the flux line (eg, |y | >> cr, or ±X > 0), as shown previously 
by Olariu & Popescu (1983, 1985). For fast wavepackets which hit the flux line, taking 
-Ptrans(^) to bc of ordcr l/cr for |y| ~ a, we get that the impulse is of order h/a. 
Therefore, it is of the same order as the dispersion in momentum, in contrast with the 
case of slow wavepackets, for which the impulse is much smaller (by a factor of a/R) 
than the dispersion. 

The ^/-component of the forward impulse, ie the imaginary part of (0), is given in 
the forward case by 

(^1 /+, \^) = ±;isin27rae(±X)Ptrans(2/). (63) 

This can be regarded as an analogue in the time domain of Shelankov's formula ( PP]) for 
the transverse momentum imparted to a stationary paraxial beam. 

5. Enclosed and distributed fluxes 

Two well-known regularizations of the Aharonov-Bohm flux line are to enclose the flux 
in an impenetrable cylindrical barrier, or to distribute the flux uniformly in a cylindrical 
tube. Here we show that the force and impulse operators in both cases approach 
the Aharonov-Bohm limit, in a sense to be explained, as the radius e of the cylinder 
approaches zero. 

In a circularly symmetric gauge, the vector potential for both models is of the form 
A'^(r) = A'^{r) 0. The kinetic momentum is given by 

MV^ = MiV: + iV;^) = ye^<^ (dr + ^ + l^^'(^)) , (64) 
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and the regularized Hamiltonian by 



The eigenfunctions of the Hamiltonian and kinetic angular momentum, with energy 
E = h^k'^ /2M and kinetic angular momentum mh, are of the form 

Xlmir) = RlmirW"''^- (66) 

The radial eigenfunctions Rl^i''^) taken to be real and normalized, like the 
Aharonov-Bohm radial eigenfunctions J\m-a\{kr), according to 

T 

These conditions determine the radial eigenfunctions up to an overall sign, which is fixed 
by requiring that Rlmi''^) approach J\jn-a\{kr) as e approaches zero. 

Let Cm{k) denote a smooth, normalized probability amplitude for energy and 
angular momentum. Let ip{r) and tp'^{r) denote the corresponding coordinate 
wavefunctions for the Aharonov-Bohm and regularized Hamiltonians, respectively. That 
is, 



R;,mir)RlJr)rdr='-^^. (67) 



^(^) =^ J2 I cUk)Jim-a\{kr)e'"''t'kdk, (68) 

m=—oo " 
1 °° poo 

^^(^) = ^ E / Uk)Rlr^ir)e'"^'f'kdk. (69) 

m=—oo ^ 

From the preceding discussion, it is clear that ip^{r) approaches ip{r) pointwise as e 
approaches zero. It turns out that their force and impulse expectation values also 
coincide as e — 0, ie 

yimmr'\r) = {i'\:F\ij), (70) 

= {m±\i^) ■ (71) 

Note that (^) and ( [7I| ) do not imply, nor is it the case, that the operators JF^ and 
approach their Aharonov-Bohm counterparts, and X±, as e approaches 0. Indeed, 
neither does the regularized Hamiltonian approach the Aharonov-Bohm Hamiltonian 
H] given e > 0, one can construct wavefunctions whose energy expectation values with 
respect to H'^ and H differ by arbitrarily large amounts. 

Instead of ( [70| ) and (|7ID, we show below, for the enclosed and distributed fluxes 
separately, that the eigenstate matrix elements of the regularized force operator 
approach the Aharonov-Bohm limit as e 0, ie, 

jim(Xp,„| \Xk,rr) = {Xp,n\^\Xk,m) = —siunak^p^ °5m,a5n,a+l- (72) 

Formally, of course, ([TOp and ([7^) are equivalent. However, for the sake of brevity 
we shall omit the details required for a rigorous demonstration. These details are 
straightforward to supply, and are similar to those given in the first part of the Appendix. 
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The result (^) for the impulse follows from the corresponding result (^) for the 
force, once it has been established that the force expectation values {ip'^l J^^{t) and 
{'4)\ J-'{t) are integrable in time. In the Appendix it is shown that, in the Aharonov- 
Bohm case, the force expectation value decays as a similar argument may be given 
for the regularized force. 



5.1. Enclosed flux. 

The kinetic momentum MV^ and Hamiltonian H*^ have the same operational form as 
in the Aharonov-Bohm case, but act on wavef unctions defined on r > e which vanish 
at r = e. For two such wavefunctions, ip''{r) and r]'^{r), assumed to be smooth and 
normalized, we have 

(^'1 J^' \ri') = 4- \MV'ri') = I UH'ip' \MV'ri') - \MV'iH'ri'))] 
at n 

= ^(^ {H'ljj' \MV'7]') - {i)^ \H' {MV'r]')) ) . (73) 

The last equality follows from the fact that the commutator [H'^,MV'^] is proportional 
to the Lorentz force operator (|I]), which vanishes for the enclosed flux. However, the 
final expression in (^) does not vanish; the relation {H'^ip'^ l^"^) = {tp'^ \H'^^'^), where 

f(r) = (MVV)W, (74) 

need not hold, because C,'^{r) need not vanish at r = e (alternatively, 1^*^) is not in the 
domain of H^). Indeed, integration by parts in ( |65D gives 

{H-r 15') - ii.- m-} = -iul I [(c + V + ' ) 



r dr ( 



2M 



'\ri^,<Pme,<P)d<P (75) 



From (H) and (^, 



r(e,0) = -e'V(6,0). (76) 
1 



Substituting this result into (|7^), we get 

m^'m = ^e^'^^r(e,0)C(e,0)e'*d0, (77) 

a result obtained previously by Peshkin (1981, 1989). Note that if we were to substitute, 
for ip'^{r) and 't]'^{r), the leading-order behaviour (|19D of Aharonov-Bohm wavefunctions, 
we would recover, formally, the Aharonov-Bohm result (|2l|) for the force expectation 
value. 

Instead, we take ip''{r) and ri^{r) in (^) to be eigenfunctions of the regularized 
Hamiltonian. Then 

(Xp,n| \ Xk,m) = ~^^-Rp,m+l'(^)-Rfc,m'(^)'^n,m+l- (78) 
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(Strictly speaking, this is not legitimate, as there would appear boundary terms at 
r = cxD in (1751). However, these would vanish when we consider expectation values, as 
in (|7g).) 

To evaluate ( [78D we need the derivatives of the radial eigenf unctions at r = e. The 
radial wavefunctions themselves are given by 

Rk,mir) = Clm {N\m-a\{ke)J\m-a\{kr) - J\m-a\{k(^) N\m-a\{kr)) , (79) 

where Ny{z) is the Neumann function. The constant is determined by the 

normalization condition (|67D , and is given by 

Clm = {J\m-a\ike) + iVf^„.| (A^e)) , (80) 
and, to leading order in e, by 

Ct^ = \Nrn-a{ke)\-' = -—^^ 4) . (81) 



r(|m-a|) V 2 

The Wronskian relation, J^{z)Nl{z) — Jl{z)N^{z) = 2/{Trz), implies that 

Substituting ( |8TD and (|8^ ) into (fTHD, we obtain, to leading order in e, 

{Xp,n\ > \Xk,^) - M r(|m + l-a|)r(|m-a|) ^^^^ 

In the limit e — 0, only the m = a term survives, and the reflection formula for the 
F-function gives 

r / *^ I I \ — ^^^^ " 

^,l^\Xp,nK " M 1(1 - a)r(a) '"''^ 

= ^sinvrafcV "(^m.a^n.a+l, (84) 

in accord with (|72p. 



5. ^. Distributed flux 

The distributed flux model was used by Nielsen and Hedegard (1995) to obtain, from 
the force balance equations, the on-shell matrix elements of the force in the limit e ^ 0. 
Here we carry out a different calculation to obtain the general matrix elements of the 
force. 

It suffices to consider the case a > (the case of negative flux is obtained from 
time-reversal). The vector potential is given by 

A'{r) = a$or/(27re2), r < e, (85) 
= a$o/2vrr, r > e, (86) 

corresponding to the magnetic field B'^{r) = (a<l'o/7re^)6(e — r), where Q{x) is the 
unit step function. In this case, the force operator is just the Lorentz force (|l]). It 
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is convenient to introduce the dimensionless radial coordinate u = r^/e^, so that the 
interior of the flux tube is given by < u < 1. The kinetic momentum is given by 



2du + '-^ + a 

u 



(87) 



Then 



T= -i 



2Mc 
2h^ 



Me3 



6(1 -m) ( 29„ + i^ + a 

u 



6{u 



The matrix elements of the force are given by 

{xl,n\ |Xfc,m> = -27r— a5„,m+i 



u 



T, 



M2 d-U - t. 



p,m+l fc,m 



(89) 



where T^^i'^) denotes the radial eigenfunction expressed in terms of the scaled variable 
u, and tl „^ = T^,^(l). 

Inside the flux tube, the radial eigenfunctions are given by (Landau & Lifschitz 
1965) 



C. 



-au/2^\ra\/2j^ 



(ke) 



+ 



\m\ — m + 1 



m\ + l,au), 0<u<l, (90) 



4a 2 

where M(a, b, z) is the confluent hypergeometric function (Abramowitz & Stegun 1970). 
Outside the flux tube, 



R 



'k,Tn \ 



Dlm-J\m-a\{kr) + El^N\m-a\{kr 



r > e. 



(91) 



The coefficients ^, Dl ^ and ^ are determined by requiring the radial 
eigenfunction and its first derivative to be continuous at r = e (the second derivative is 
then continuous there as well, as it turns out), and by the normalization condition 

{DIJ' + {EIJ' = 1, (92) 

which follows from (p7|). 

To evaluate the force matrix element (|89D, we only require the function inside the 
fiux cylinder. Straightforward algebra gives the coefficient C^^^, to leading order in e, 

as 

2e°/2(lA;e)l'"-"l 



r(|m - a|) [{\m - a\ + \m\ - a)fm + 2/4 



where 



(93) 
(94) 



Fm{u) = M(i(|m| — m + 1), \m\ + 1, an), 

and fm and denote the values of Fm and F^ at u = 1. 

we find that that {Xp.n \ |Xfc,m) is of order 

a. We obtain 

(95) 



Substituting (pO| ) and (0) into 
g|m-Q|+|m+i-Q|-i^ therefore vanishes in the limit e ^ unless m 

2^2 



M ^ R{a 



ni,a'-'n,a+l ) 
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where 

L{a) = -2a F'^F^+^e'^'^^-^^^+Uu + afj^+i, 

R{a) = 2{p{l)U, + fU,fa), (96) 

and 

p{u) = a + 1 — au. (97) 
As we show below, L{a) = R{a), or, equivalently, 

/' 2aF^F,+ie"(i-")n'^+i du = afja+i - 2(p(l)/.+i + f,^,)C (98) 

With this identity, (|95|) gives the required result (|72D . 

To establish the identity (|98D , it is convenient to express Fa+i in terms of Fa by 
means of the recurrence relation (Abramowitz & Stegun 1970) 

(a + |)M(i, a + 2, m) = (a + 1) (M(i, a + 1, m) - M'(i, a + 1, m)) , (99) 
which implies that 

{a+\)Fa+i = {a + l){Fa-F'Ja). (100) 
With the differential equation 

uF': = -p{u)F'^ + (101) 



W{u) = 2^«'^+ie"(^-") - p{u)FaF^ - uF'^) . (102) 



it is straightforward to show that the integrand on the left-hand side of ( pH]) is given by 
W'(u), where 

,a + 1 

W{0) vanishes, whereas Vr(l), with the aid of ( p.00| ) and (|101|) , is seen to be equal to 
the right-hand side of (pH]). 
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Appendix A. Wavepacket expectation values 

The force and impulse due to an Aharonov-Bohm flux line can be calculated rigorously 
for suitably well-behaved wavefunctions ipii^)- We will take these to be such that 

Cm{k) = {xk,m li^) is smooth in k and falls off, along with its derivatives, 

faster than any power of k and m. 



Using standard arguments, one can show that (|A.1|) implies the following properties of 
?/'(r) and {Hilj){r), where H is the Aharonov-Bohm Hamiltonian: 

ip{r) and {Hil)){r) are smooth for r > and fall off, along with their (p^2) 

derivatives, faster than any power of r. 
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and 

{Hij)Ur) = D^r\^~^\ + O (rl™--l+i) , (A.3) 
where, in general, 



^m{r) = —j^ r7(r,0)e-™*d0. (A.4) 

(In fact, properties (|A.2|) and ( |A.3| ) are also shared by {H^ilj){r), for j > 1. The 
argument to follow would hold under weaker conditions, but we assume ( [A.l| ) for 
simplicity.) 

The expectation value of force is given by 



^ /" /" - {HiP*){MVip) +ip*{MVHxlj)^ dV, (A.5) 



where MV is given by (|1^). From ( |A.2| ) and ( |A.3| ), it is evident that the r-integral 
in ( [A.5| ) converges absolutely. This allows us to introduce a factor exp(— e^r^) in the 
integrand, and then take the limit of the integral as e 0. This Gaussian factor will 
justify subsequent reorderings of operations. Note that the integral cannot be expressed 
in terms of the expectation value of the commutator [H, MV], because of the singularity 
in the radial derivative of iplr) at the origin (specifically, {MVip){r) is not in the domain 
of H). 

We introduce the eigenfunction expansion 

1 °" poo 

V'(r-) = — Yl / Cmik)Jim-a\ikr)e'"''^kdk, (A.6) 



m=— oo 



and a similar expansion for {Hip){r), with Cm{k) replaced by —{fi^k^ /2M)cm{k). Using 
standard arguments, one can show that (|A.1|) implies that the differential operator 
MV, when applied to ip and Hip, can be taken inside the m-sum and fc-integral. The 
recurrence relation, 

Ju±i{z) = t[jI{z)t^M^)), (A.7) 

implies that 

MV {J\^_^\{kr)e''^^) = sgn{m - a)ihkJ\m+i-a\{kr)e'^"'+^^'^ , m ^ a, 

= -ihkJa-i{kr)e'^''+^^'^, m = a. (A.8) 

Substituting ( |A.6D and ( |A.8| ) into (|A.5|) , along with the eigenfunction expansion of 'ip*{r) 
with coefficients c*(p), we obtain 

fc2 /"OO p2tt °° °° POO POO 

— / e-^'"Vdr / d0 V S^Qiim+i-n)^ dp / ^j. 



m=— oo n=— oo 



xc:{p)Uk)k'p{e-p')J\r.^^\{pr)x{ ^/ / l'-^' ^' ^ ' (A.9) 

I —Ja~i{kr), m = a. 
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The sums and integrals in ( |A.9| ) are uniformly and absolutely convergent, and can 
be interchanged. On performing the 0-integral, the sum on n collapses to the single 
term n = m + 1. We obtain 

ft" 



A-kM 

where, for m ^ a, 

K^{k,p) = sgn(m — a)/c^p(/c^ — 

and, for m = a, 

Kl{k,p) = ep{p'-k') [ 

Jo 



roo POO 

hm ^ / / Kl{k,p)cl,^Mcm{k)dMp, (A.IO) 
Jo Jo 

m=— CO " " 



' Jy{pr)Jy{kr)r dr, z/ = |m — a + 1|, (A. 11) 



Ju{pr)J-u{kr)r dr, z/ = 1 — ct. (A. 12) 



Below, in Appendix A.l , we show that the contributions from the m ^ a terms vanish 
in the limit, ie 



lim > 

e^O ^ 



OO /"OO 



Jo 



KUk,p)c*{p)Ca{k) dkdp = 0, 



(A.13) 



while, in [Appendix A.2|, we show for the m = a term that 



lim 



OO poo 



JO 



Kl{k,p)cl^Mca{k)dkdp 



c2 roo poo 

— sin Tca 

-/O JO 



A;i+V-"d(p)c„(A;)dMp. 



(A. 14) 



Substitution of i\KA^ and (|XA^ ) into ([AA0|) gives 



T\ 



smvra 



OO roo 



^0 



k^^y'-'^cl^Mca{k)dkdp. (A.15) 



This is equivalent to the matrix element (p5|), obtained formally in Section ^ 

We note that, with e = 0, the integrals in ( |A.ll ) and ( A.12| ) correspond to singular 
(ie, not absolutely convergent) cases of the discontinuous Weber-Schafheitlin integral 
(Abramowitz & Stegun 1970). Formal evaluation of these integrals would give ( |A.13| ) 
and ( |A.14|) immediately. The arguments in Appendix A.l and Appendix A. 2 serve to 
justify these formal results. 

To obtain the force expectation value (plD, we express Caik) and c*_^]^(p) in ( |A.15| ) 
in terms of ip{r) to obtain 

2h^ 



sin Txa lim 

M 



^l+ap2-<5g-e2{fc2+p2) 



^0 

X ( / i^*a+l{s)Jl-&{ps)sds 



il)a{r)Ja{kr)r dr ) dkdp 



(A.16) 



Note that the convergence factor exp(— e^(fc^ + P^)) can be introduced, and the limit 
e — s> taken outside the integral, since, by ( |A.1D , the k- and p-integrals in (|A.15|) are 
absolutely convergent. By ( |A.2|) and (|A.3|) , the r- and s-integrals in ( |A.16| ) are absolutely 
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convergent, so that, for e > 0, we can interchange the order of integration. The k- and 
p-integrals can be evaluated using (pHf), with the result 



M e^o Jq \ J v4eV \4eV 



X 







smTr&CaC*^iT{l + a)r(2 - a) 



M 

= - (A.17) 

where the coefficients Ca and Ca+i are given in ( |A.3| ). This is just the result (|2T| ) of 
Section 0. 

Concerning the impulse, it is straightforward to justify, using arguments like the 
preceding ones, the calculations of Section |] leading to (^61). It is only necessary 
to check that the time-dependent expectation value, {tp(t)\J-'\ip(t)), is integrable in 
t. {ip{t)\ \ip(t)) is given by an expression like ( [A. 151 ), but with Ca{k) and c*^i(p) 
modulated by the factors exp{—ihk'^t /2M) and exp(i^p^t/2M) respectively. We have 
that 

j%^+^c^{k)ex^{-\^t^ Ak 

s(*^^».W)exp(-i^*)d. = 0(lA) (A.18) 

for large t (the integration by parts is justified by ( |A.2| ) and ( [A. 3D ). Similarly, 

p'^-'^cl^^ip) exp{ihpH/2M) dp = 0(l/t). It follows that {ip{t)\ \ilj{t)) falls off as 

Appendix A.l. Proof of l\A.l'^ 

Given functions / and g defined on a domain D, we will say that / is dominated by g 
if, for some constant C, |/| < Cg throughout D. For functions indexed by an integer 
m, eg fm and g^, we will say that fm is dominated by g^ if \fm\ < Cgm for some C 
which does not depend on m. Thus, from (|A.1| ), 

Cfn{k)c*^+i{p) is dominated by (1 + m^)~^{l + k'^ + p^)'^ for k,p > 0. (A.19) 

The integral in ( |A.11| ) can be evaluated (Gradshteyn and Ryzhik 1980) to give 

where Iu{z) is a modified Bessel function, and u = \m — a\. From the asymptotic 
behaviour of Iu{z) for large argument, it follows that Iu{z) is dominated by e^/A/i for z 
real. Therefore, the left-hand side of (|A.13|) is dominated by 



Force and impulse from an Aharonov-Bohm flux line 



22 



Let us divide the domain of the (A;,p)-integral into regions inside and outside the strip 
|fc — p| < e^/^. Inside, the integrand is dominated by e'^^^k^/{l + k'^)^; thus the integral 
over the strip is dominated by e^/^. Outside the strip, the integrand is dominated by 
exp (— e"^''^) /(I + /c^ thus the -integral outside the strip is dominated by 

exp (— e"^'''^) . Therefore, the expression in ( |A.21| ) is dominated by 



^ 1 1 



^-^ 1 + e 

m=—oo 

which vanishes as e ^ 0. 



Appendix A. 2. Proof of l \A.14 ) 



Substituting the series expansion 



2/ ^ uWiu + v+l) \ 4 

n=0 ^ 



and a similar expansion for J-u{z), we get that 

POO ^ ^ 

I e~^ Ji,{pr)J^^{kr)r dr 



S^vv nu+v+i) ( k^vf P^y 

For e > 0, the r-integral and u— and f-sums are absolutely convergent for all k,p > 0. 
Inserting in ( |A.24| ) the integral representation for the reciprocal of the beta-function 
(Gradshteyn and Ryzhik 1980) 

r{u + v + i) _ 1 

r{u-u + l)T{v + U + ~ {u + v + l)B{u -u + l,v + u + l) 

= -Re / (2isinre-'")" (-2isinre'")%2^''("-"/2)dr, (A.25) 
we may perform the sums to obtain 

POO 

k" / e"" ^ Ju{pr)J^u{kr)r dr 
Jo 

= -^Re exp - -^ , sin^r + i^ ^ sin 2r + 2iz/(r - 7r/2) dr. (A.26) 

Substituting this result into ( |A.12| ), we get 

Kl{k,p) = — A;1+V""(P^ - fc^)Re / e"^/^' dr, (A.27) 
T^e Jo 

where the exponent S is given by 

^ = 1 +p2) (1 _ cos 2r) - ii(p2 _ sin 2r - 216^ (1 - a) (r - n/2) . (A.28) 

It is clear that the main contribution to the r-integral in ( [A.27] ) comes from the 
neighbourhood of r = 0. If S* is expanded about r = to second order, the r-integral 
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yields an error function, whose leading-order asymptotics as e ^ leads directly to 
the required result ( [A.14| ). However, the next term in the asymptotic expansion is not 
uniformly bounded in k and p - it contains a factor {k"^ —p^)""^ ^ so we must take some 
additional care. 

To proceed, we divide the domain of the A;,p-integral into the three regions specified 
below, writing the left hand side of ( |A.14 ) as 



Ihn j +lim j +lhn j j Kl{k,p)cl^,{p)ca{k) dkdp, (A.29) 



and analyzing the contribution from each region separately. 

Let Di denote the region k,p > 0, k"^ + p"^ < e^, where /3 is chosen to satisfy 



I<f3< 



3' 



(A.30) 



In this region, the exponential factor in ( |A.27| ) is bounded, and k^^"'p'^~°'{p'^ — k"^) is 
dominated by e^^/^. The coefficients |ca(fc)| and |c*^i(p)| are bounded, and the area of 
Di is dominated by e^^, so that the contribution from Di in ( [A.29|) is dominated by 
g7/3/2-2^ Given the choice of /?, this vanishes as e ^ 0. 

Let D2 be the region k,p > 0, k"^ + p"^ > and {p"^ — k'^\ < e'', where 7 is chosen to 
satisfy 



(A.31) 



since /? < |, the inequality (|A.31|) can be always be satisfied). Since 1 — cos2r > for 
< r < 7r/2, the factor exp(— S'/e^) is dominated by the Gaussian exp(— cr^r^), where 

Then, from i\K27l) , 



a 



1^/3/2-1 



Thus JJ^^^ exp(— S/e^) dr is dominated by e 



1-13/2 



Kl{k,p) is dominated by e7-/3/2-i^i+«p2-a jj^ p^^^ ( |OD , p'^-°'cl^^{p)caik) is 



integrable over the region k,p > 0. Therefore, the contribution from D2 to (|A.29| ) 



is dominated by e '^'2'^^ (the additional factor of e'~' is due to the fact that the integral 
is confined to \k'^ — p'^\ < e^). Given the choice of 7, this vanishes as e ^ 0. 

The remaining region D3 is given hy k,p > 0, k^ + p"^ > and Ip^ 
Integrating by parts with respect to r in (|A.27|) , we get that 



k^\ > e\ 



K{k,p) = -k 

n 



l+y2-a^^2 _ j^2 



X Re 











y Sr 


r=0 


Sr 


-j 

r=iT/ Z 



f7r/2 



: Sr- 

~sl 



■dr 



(A.32) 



The first term gives 
2 

TT 



-sin7raA;^+V~" 



1 + 



4(1 



«)e2 



(A.33) 
0, the required 



p2 _ J^2 

Its contribution to the integral over in (|A.29| ) yields, in the limit e 
result (prT^ ). 

It remains to show that the contribution to the Ds-integral from the remaining 
terms in (|A.32|) vanishes in the limit e — 0, It is readily seen that the contribution from 
the second term vanishes exponentially with e. For the third term, we note that, on 
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the interval < r < 7r/2, S^r/S'^ is dominated by +p^)/(/c^ —V^Y-, and exp(— S'/e^) 

is dominated by exp(— e'^^^r^/2). Therefore, the integral JJ^''^ e~'^/^^(S'T-T-/S'^)dr is 

dominated by e^~^/^(A;^ ~ A;^)^. Thus, the third term in (|A.32|) is dominated 

by e^-^/2(^A;2 +j92)^i+«p2-fi/|p _p2|^ ^j^j^j^ dominated by e^~^/2~'^(A;2 +p2)7/2_ 

The contribution to the integral over D3 in (|A.29| ) is dominated by g^'/^/^-T^ which, by 
the choice of 7, vanishes as e — >^ 0. 
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